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Strong isospin breaking in the spectrum of the nucleons and deltas can be studied in lattice QCD 
with the help of chiral perturbation theory. At leading order in the chiral expansion, the mass 
splittings between the proton and neutron and between the deltas are linear in the quark mass 
difference. The next-to-leading order contributions to these splittings vanish even away from the 
strong-isospin limit. Therefore, any non-linear quark mass dependence of these mass splittings is a 
signal of the next-to-next-to-leading order mass contributions, thus providing access to LECs at this 
order. We determine the mass splittings of the nucleons and deltas in two-flavor, heavy baryon chiral 
- - - perturbation theory to next-to-next-to-leading order. We also derive expressions for the nucleon and 

delta masses in partially quenched chiral perturbation theory to the same order. The resulting mass 

. expressions will be useful both for the extrapolation of lattice data on baryon masses, and for the 

' study of strong isospin breaking. 

{3JQ' PACS numbers: 12.38.Gc, 12.39.Fe 
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I. INTRODUCTION 



Understanding nuclear physics from first principles is a long standing challenge. Quantum chromodynamics (QCD) 
is the well established theory of strong interactions that describes the underlying interactions of colored quarks and 
gluons which bind together to form the color neutral hadrons observed in nature: nucleons, pions, etc. The equations 
of motion for the quarks and gluons are non-linear, and at energy scales relevant for nuclear physics, _E ~ 1 GeV, the 
theory is non-perturbative. Thus the description of hadrons from first principles remains elusive. One can write down 
an effective theory of QCD, chiral perturbation theory (^PT) 1], which provides a model-independent description 
of the low energy interactions of the hadrons by preserving the symmetries of QCD. xPT is written in terms of the 
pions (pseudo-Goldstone bosons which appear from spontaneous chiral symmetry breaking) and their interactions 
with baryonic states, like the nucleons and delta-resonances. As the pions are the lightest particles in the theory, 
they dominate the low-energy interactions of the hadrons. In xPT, short distance physics is encoded in an infinite 
'] tower of higher dimensional operators, whose coefficients are low-energy constants (LECs). The utility of xPT is that 
t these higher dimensional operators are suppressed by appropriate powers of the cutoff, where operators with higher 
^ inverse powers of the cutoff are higher order in the expansion. Thus, to a given order, one needs to keep only a finite 
number of operators, and the theoretical error associated with a given calculation is indicated by the size of omitted 
terms. Symmetry constrains the number of operators which enter at a given order, but not the numerical values of 
their coefficients. Thus the LECs must be fit by comparing either to experiment or lattice QCD. 

Lattice QCD is a tool one can use to numerically calculate QCD observables from first principles — the mass 
spectrum of baryons and mesons being notable examples. Due to computational restrictions, however, the light 
quarks must presently be simulated with unphysically large masses. So that observables calculated on the lattice 
can be extrapolated from the larger quark masses used in lattice simulations to the the physical quark masses, one 
requires an understanding of the quark mass dependence of QCD. For low-energy QCD, xPT is precisely what is 
needed, as the quark masses explicitly break chiral symmetry. For the up and down quarks, this breaking can be 
treated perturbatively; thus, these quark masses enter xPT as parameters of the theory, allowing QCD observables 
to be calculated as a systematic expansion in powers of the quark masses. Consequently lattice QCD can be used in 
conjunction with xPT to first determine the LECs relevant to a given observable, and once these LECs are determined, 
extrapolate to make rigorous first principles predictions of QCD observables. 

Lattice QCD calculations of the baryon masses began in the early 1980's but until recently, unquenched lattice 
calculations have been limited to quark masses such that the pion has been rather heavy, to ^ > 500 MeV ^ ^, ^, ^ 
III 1^ l2l lldl • There are a handful of fits to the lattice data on baryon masses using xPT [nllialiauj. However, the 
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viability of xPT with such large pion masses is questionable |l5lll6| . Ideally, one would like to see lattice data with all 
the pion masses less than ^ 400 MeV, to properly map out the chiral regime. This is a very exciting time, as we are 
beginning to see lattice data of baryonic observables with TOtt ~ 330 MeV, e.g. 0. If not the pions accessible today, 
then the ones in the near future will be light enough that one can apply chiral perturbation theory to observables 
calculated from lattice QCD. We anticipate that lattice QCD in conjunction withyPT and PQxPT will soon be 
making predictions of nuclcon observables with theoretical errors at the 20% level llfil.^ 

Calculation of the baryonic mass spectrum from lattice QCD data in the chiral regime will be a first step in 
understanding nuclear physics from first principles. To this end, we determine the masses of the nucleons and 
delta-resonances to next-to-next-to-leading order (NNLO) in two-flavor non-degenerate xPT. We highlight the mass 
splittings among the nucleons and deltas, as they provide a way to determine some of the LECs entering at NNLO. 
This is because strong isospin breaking in the baryon spectrum receives the first non-linear quark mass dependence 
at NNLO. In the appendix we calculate the masses of the nucleons and deltas in PQxPT, with two non-degenerate 
valence and sea quarks. These expressions are necessary to extrapolate partially quenched lattice data of the nucleon 
and delta masses to the physical values of the quark masses. 



II. MASSES IN xPT 

In this section, we calculate the masses of the nucleons and deltas in the non-degenerate SU(2) flavor group. We 
start by reviewing the chiral Lagrangian for mesons and baryons. Next we calculate the masses of the nucleons and 
finally the masses of the deltas. The masses of the nucleons and deltas have been investigated considerably in xPT, 
see e.g. lUi, ilS, i21i i2&, 111] • In Appendix 1X1 we take up the partially quenched calculation. 



A. xPT Lagrangian 

For two flavors of massless quarks, the QCD Lagrangian exhibits a chiral symmetry, SU(2)L®SU(2)i?(g)U(l)y. This 
symmetry is spontaneously broken by the vacuum down to SU(2)v(E'U(l)y . xPT is the effective theory generated 
by expanding about the physical vacuum state of QCD. Without the source of explicit chiral symmetry breaking 
generated by the quark mass term, the pions would appear as the Goldstone bosons of this symmetry breaking. In 
reality the up and down quark masses are not zero but are quite small compared to the scale of chiral symmetry 
breaking. Thus the pions appear in nature as pseudo- Goldstone bosons. 

To build the effective theory, we collect the pseudo- Goldstone bosons in an exponential matrix 

2 i <I>\ _2 




S = exp(^-^j =e^ $-11^- 1^0 I- (1) 

The pion decay constant / is 132 MeV in the above conventions. The effective Lagrangian describing the dynamics 
of the pions at leading order in xPT is 

ytr(a''E^9^E) +Atr(77iJS + m,E^) . (2) 
Above the quark mass matrix mg is given by 

m, = diag(m„,md). (3) 

Expanding the Lagrangian Eq. (01 one finds that to leading order the pions are canonically normalized with their 
masses given by 

2 4A, , 
^ j2{^u + md). (4) 



^ We mean that the xPT and PQxPT expressions will be fit only to lattice calculations but not experimental values. These expressions 
will then be used to extrapolate down to the physical quark masses, make predictions of experimental observables, and ultimately 
predictions of observables not currently accessible via experiment. 
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The baryons are three quark states in the spectrum of QCD. To include the lowest-lying spin-i and spin-| baryons 
into xPT, we use heavy baryon xPT (HBxPT) 0, yO, |^ 113, 113 . As is well known, this approach leads to a 
consistent power-counting scheme for including the baryons. In SU(2), there is a doublet N of spin-i nucleons 

»=(f.). (^) 

and a quartet of spin-| deltas. These latter states are contained in the completely symmetric tensor T*^*^, where 

r"i = A++, r"2= 1a+, ri22^^A", and t222 = a-. (6) 

•\/3 V 3 

The free Lagrangian for nucleons and deltas to leading order in HB^PT can be written as 

£ = iNv-DN + 2aMNMN + 2aMNNtr{M) 

-if'' vDT^ + AT^T^ + 2^mT^ MT^- 2aMT^T^ tr(X). (7) 

Above, A is the leading-order mass splitting between the nucleons and deltas in the chiral limit. As a dimensionful 
parameter it must be treated in the power counting and we treat A ^ 0{q) [34 .35.] . The chirally covariant derivative 
acts on the nuclcon field as 



iD^N),^d^N, + {V''y N,,, (8) 

and on the delta field as 

{D^T),jk = d^T,,k + {Vn/T,,,k + {Vn/T,j,k + Cl^^)fe"'ryfe' • (9) 
The vector and axial-vector meson fields are given by 



and finally the spurion mass field 



M^^i^^m,^^ +^ml() . (11) 



The Lagrange density describing the interactions of the nucleons and deltas with the pions at leading order can be 
written as 

£ = 2gANS ■ AN + gAN {T^^A^N + NA^T^) + 2gAAT'' S ■ AT^, (12) 

with as the covariant spin operator. Additionally we need higher-order terms to calculate masses at NNLO. The 
first set of higher-order terms are all operators that are constrained by reparameterization invariance 

These 

terms must be included to insure the Lorentz invariance of HBxPT at 0{1/MbY , where Mb ~ A,^ is the average 
nucleon mass in the chiral limit. As a result of reparameterization invariance, these operators have fixed coefficients 
and explicitly these operators are 



^ For an earlier, alternate approach for the nucleons, see, e.g., Ref. |32| . 
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where Dl = - {v ■ Df. 

Next there are three sets of higher-dimensional operators whose coefficients are not pre-detcrmincd. The first set 
makes contributions to the mass of the nucleons at NNLO. The operators in this set are 

£ = -^Ib'l'N M\N + hl'l^NiY{M\) + hl'N MNtY{M) + b^'NN[tY{M)f 



-b^NNtiiA ■ A) + b^^NNtriv ■ Av ■ A) 



(14) 



The LECs bf^, b^, and b""^ are aU dimensionless. The choice in numbering the coefficients was made to be consistent 
with as wiU become clear in Appendix^ where we address the partially quenched calculation. The second set of 
higher-dimensional operators with undetermined coefficients consists of operators that contribute to the delta masses 
at NNLO. These operators are 

^ = ^[t2TfiA^^V'iAn/Tl:^,,+4 {T,T^)tr{A,An 

+t§TfiAn/{A,)/TI1^,, + 4T,n{A>^A,) 
+q^Tf{v ■ A)/{v ■ A)/tP^,, + e/T,T^Mv -AvA) 

+tf (T^XT'^) tr(X)+ if T^r^[tr(X)]2l (15) 



and all of the LECs if^, tf-, t"^^, and tf are dimensionless. 

The last set of higher-dimensional operators with undetermined coefficients all involve the mass-splitting parameter 
A, which is a singlet under chiral transformations. Including the spin-| fields in xPT thus requires the addition of 
operators involving powers of A/A^. We shall not write out such operators explicitly. To account for the effects of 
these operators, all LECs in the calculation must be treated as arbitrary polynomial functions of A/A^ and expanded 
out to the required order. For example 

/A\ / A A2 

VAx/ V A;, A2 

A \ A 



Determination of these LECs requires the ability to tune the parameter A and for this reason we do not keep such 
operators explicitly. Furthermore at this order the nucleon mass, Mq, and the nucleon-delta mass splitting in the 
chiral limit are also functions of A/At^. This variation, however, can be absorbed into a redefinition of the parameters 
Mo, and A (including contributions from loop graphs). We use these value for Mq and A everywhere. Differences 
resulting from using this value of A in expressions for loop graphs are beyond the order we work. 



B. Nucleons 



The mass of the i*'' nucleon in the chiral expansion can be written as 

Mb, = Mo (A) - M^ (^, A) - M^ /'^ (^, A) - M^f (a*, A) + . . . (17) 

Here, /i is the renormalization scale and A is the quark mass independent delta-nucleon mass splitting. Mq (A) is the 
renormalized nucleon mass in the chiral limit, is independent of fi, of and also of the Bi. Afjj"^ is the contribution 
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to the i*'' nucleon of the order m 



(") 3 



Throughout this work, we use dimensional regularization with a modified 



minimal subtraction (MS) scheme, where we have consistently subtracted terms proportional to 

- - 7b + 1 + log47r. 

e 

The diagrams relevant to calculate the nucleon mass to NNLO are depicted in (2^. We find that the leading-order 
contributions to the nucleon masses are 



Mg^ = 2aMmB + 2aM ti{mq), 



(18) 



where 



ruB 



ma, B ^ n 



(19) 



89. 



The next-to-leading contributions are 

,,(3/2) _ 3 2 3 8g|jv 

- WttP^^"^- + 3(4^/)2-^^"'^'^'^^' 
where the function J- is defined by 

!F{m,S,ii) = (m^ — (5^) 



6 — -i/ 6^ — im? + ie 



\/ — vn? log 



(51oe 



Finally the next-to-next-to-leading contributions to the nucleon mass are 



M 



(2 



1 



g - [Zb- l)M)^' + — {6f {ruB? + 6f tr(™2) + foM ,^.(^^) ^ [tr(™J]2} 



1 



(47r/)2 
36^ ^ 



Cf/:(m,,Ai) 



(47r/)= 



(4^/)2 
3 6^^ 



4(4^/)3 



+ 



1 



(47r/)2 16AfB 



(4^/)2 

(47r/) 



tr(TOq) 



6 
2 



1 55 



2 



(47r/)2 2Mb 



(4^/)2 



tr(mg) [j'(m^. A, ^) + ml] 



2 ^ 7r 



2 

£{m^,n) + -m^ 



AW 



jmG^ [j'(m^, A,/i) 



(20) 



(21) 



(22) 



Here we have used (mB)^ which is the square of the tree-level coefficients appearing in Eq. p9|l . Above the 
wavefunction renormalization Zb is given by 



9g 



2(47r/)2 



C{m^,fi) + -ml 



(4^/)2 



[Jim,,A,fi) 



(23) 



The renormalization scale dependence appears at each order because we are treating the LECs as polynomial functions of A. If we 
instead treat this dependence explicitly and expand the nucleon mass in powers of q, where A ~ 5 and mg ~ q^, the mass then takes 
the form 

»,r (2) (3) (4) , 

Mg^ = mg — rrig' — rrig' — mg' + . . . , 

with rrSg^ as the renormalization scale independent mass contribution to the j"' nucleon strictly of the order 5" . The parameter mp is 
nucleon mass in the chiral limit, but only to leading order in A, i.e. mg = A/o(A = 0). 
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TABLE L The coefficients , , and in xPT. Coefficients are hsted for the nucleons. 









r<B 


p 


2aM(2mu + ma) 


aM{mu + 2md) 


|(7m„ + 2md) 


n 


2aM{mu + 2md) 


aM(2m„ + vrid) 


|(2m„ + 7md) 



TABLE IL The tree-level coefficients in xPT. The coefficients ttit, (m'^)T, and (mm')T are listed for the deltas T . 





TTlT 




{mm')T 


A++ 




3ml 


3ml 


A+ 


2mu + md 


2ml + ml 


ml + 2m„md 


A« 


niu + 2md 


■ml + 2ml 


2mumd + ml 


A" 


3md 


3ml 


3ml 



The coefficients in the NNLO contribution, namely , , and G^, are given in Table |l] and depend on whether 
B = p OT B = n. The equations above, Eq. H22() and Eq. H2,S|1 . also employ abbreviations for non-analytic functions 
arising from loop contributions. These functions are 

9 

C{m,fj,) ^m^ log— , (24) 

777,^ 

>C(m,/i) = log— , (25) 



:7(m, 5, /i) = (m^ - 2^2) i^g + 2<5 V^^^ log ( ^ "^f' + ) + 2^5^ log (^] . (26) 

The expressions we have derived in this section, as well as those throughout this work, are functions of the quark 
masses; e.g. TOtt above is merely a replacement for the combination of quark masses given in Eq. Q. These expressions, 
thus require the lattice practitioner to determine the quark masses. In SU{2) xPT away from the isospin limit, there 
are no equivalent expressions in terms of the meson masses, as one cannot independently associate m„ and to^ to the 
pion masses, which are degenerate to the order we are working. One can only equate the average value, ^(m„ -I- m^), 
with the pion mass. Therefore one cannot plot the baryon masses as a function of the physical meson masses, unless 
one works in the isospin limit, m„ = m^. This is unlike the case in the isospin limit of ^[/(S), where there are only 
two independent quark masses, but three independent meson masses, and one can always convert from a quark mass 
expansion to a meson mass expansion via the Gell-Mann-Okubo relation.^ Lastly we remark that if one is using xPT 
to determine the quark masses from meson masses that are determined on the lattice, one must use the one-loop 
expression in xPT, else one looses contributions to baryon masses that are of NNLO. Thus the quark mass expansion 
is the one to use in SU{2) xPT, at least away from the isospin limit. 



C. Deltas 

The mass of the z*'' delta in the chiral expansion can be written as 

Mt. = Mo (A) + A + A4f (m, A) + M^f (p, A) + (/i, A) + . . . (27) 

Here, Mq (A) is the renormalized nucleon mass in the chiral limit from Eq. (|17|) . and A is the rcnormalized nucleon- 
delta mass splitting in the chiral limit. Both of these quantities are independent of ruq and also of the T^. Af^"-* is 
the contribution to the i*'' decuplet baryon of the order mg"\ and /i is the renormalization scale. ^ 



* This problem can be avoided in SU{4\2) PQxPT, as there are more independent meson masses than independent quark masses, so 
one can algebraically convert from the quark mass expansion to the meson mass expansion. For example, to leading order, mu = 
|i (^ml - m% + m%) . This would require one to know the mass of the mesons made of one valence and one sea quark. But if one is 
interested in the non-isospin SU{2) limit of the PQxPT expressions, the problem is unavoidable. 

As with the nucleon masses, the renormalization scale appears in each term contributing to the delta masses because we implicitly treat 
the LECs as polynomial functions of the mass parameter A. If we expand out the LECs, the /^-dependence disappears at each order. 
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TABLE III: The coefficients , , and in xPT. Coefficients are listed for the delta states T. 





Ctt 


7r 


gJ 


A++ 


2mu + rud 






A+ 


f m„ + |md 


+ ylwi'd 


^aM{2mu + rUd) 


A« 




+ 


|QAf (m„ + 2md) 


A- 




|m„ + 





The diagrams relevant to calculate the delta masses to NNLO are depicted in |27i. We find that the leading-order 
contributions to the delta masses are 



(1) 



T - -^IM m-T - 2a M tr(TOg), 



(28) 



where the tree-level coefficients niT are given in Table ITl) for the deltas T. The next-to-leading order contributions 
are 



(3/2) 



^^^^ ^(m.,-A,M). 



4327r/2 3(47r/)2 
Finally the next-to-next-to-leading contributions to the delta masses are 



(29) 



+ 4^ ("^')^ + (mm')T + tr(m2) + tut tr(m,) + tf [tr(™,) 



rC^ £(m^,/i) -f 



(47r/)2 



6 gM 

(47r/)2 



tr(mg) £(TO7r,M) + 



1 



1 



4(47r/)3 



19 4 



(47r/)3 V2 



1 



48(47r/)2AfB 



5g] 



AN 



8(47r/)2M£ 



9(4^/)2 

10gAA7A/ pT 

■ 9(47r/)2 - 



tr(mQ) 



15 



, 26 2 
£(m^,^) + — 



(^^tr(mQ)^(m.,-A,M) 



2(4^/)^ 



rG^ J^(m^,-A,/i). 



(30) 



Here we have used {m?)T, {ram')T to label the tree-level coefficients that appear in Table ITU Above the wavefunction 
renormalization Zt is given by 



Zt - 1 



25.9 



2 

AA 



18(47r/)2 



N 26 2' 



9an 



(31) 



The coefficients in the NNLO contribution, namely Cj, , and Gj, are given in Table Hill and depend on the delta 
state T. 



III. MASS SPLITTINGS 



Having derived the nucleon and delta masses to NNLO in the chiral expansion, we now focus on the mass splittings 
between these states. To begin, we consider the nucleon mass splitting, which to our knowledge was first theoretically 
addressed in js^- The degeneracy between the proton and neutron is broken by leading-order effects in the chiral 
theory, see Eq. ((TH|l . Beyond this order pion loops contribute, but to the order we are working, all the pions are 
degenerate, even away from the isospin limit. Thus the NLO contributions to the neutron and proton masses are the 
same, see Eq. ()2(J|I . and the mass splitting, M„ — Mp, is given to NLO accuracy, entirely by the difference of the LO 
mass contribution in Eq. (|18|l . and is linear in rud — Any deviation from this linear mass splitting seen in lattice 
simulations of the nucleon masses should be a signature of the NNLO mass contributions and certain LECs that arise 
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TABLE IV: The coefficient 5t , which encodes the violation of the delta equal spacing rule. 



SMt 


St'' 


M^- - Mao 


2md 


Mao - Ma+ 


mu + ma 


Ma+ - Ma++ 


2m„ 



at this order. Additionally, the nucleon mass splitting can be enhanced from that in nature by increasing the quark 
mass splitting on the lattice, nid — rriu- This enhancement, combined with the vanishing of the NLO contribution 
to the mass splitting, provides us with a means of cleanly determining the NNLO nucleon mass contributions and 
isolating certain LECs arising at this order. These effects are normally obscured by the NLO contributions. 
We find the nucleon mass splitting is given to NNLO by 



Mn - Mp 



-2aM{md - mu) 



2 

" 9 AN 



(4^/)2 

JiniT,, A, ^) 
(4^/)2 



:1M 



2aM{6gl + l) 



(6f + C) 



31 



(32) 



An identical situation arises for the deltas, as their degeneracy is broken at leading order in the chiral expansion, 
while the next contribution to their splittings occurs at NNLO. Their mass splittings are given by (5Mt, which stands 
for Ma- - Mao , Mao ~ and Ma+ - Ma++. We find. 



5Mj 



-,lM{md - m„) 



(4^/)2 
£{m„,fi) 207m 



^(^1 +^4 )-^5aa7m 



(47r/)2 



81 



27 2 



AN 



12TTf 



(4^/)2 



{iM + a A/) 



(33) 



where 5t^^ is given in Table Hvl Considering these splittings, there are a few things to note. We know from experiment 
that the neutron is more massive than the proton, and we expect that this is mostly due to rud > ruu- Similarly, 
we expect the deltas to follow a similar pattern with A~ being more massive than A'^, and so on, although this is 
presently undetermined experimentally. From these expectations and the expressions above, it is expected that aM is 
negative and jm is positive. Also, the delta masses, to a good approximation are expected to follow an equal spacing 
rule. In HB^PT, this rule is first violated at NNLO, and even then, only by the operator in Eq. I|15|) with coefficient 
. This coefficient can be isolated by taking successive differences. 

Lastly we should comment that in nature isospin violation in the baryon masses has another source of the same 
size as the NNLO chiral effects, namely electromagnetic contributions. In lattice QCD calculations one can turn off 
the electric charges of the quarks. This is the scenario for which our calculations are app licable. There have been a 
few lattice computations of electromagnetic contributions to hadronic masses |33, |41| • For a recent discussion of 
the electromagnetic effects in hadrons, see Ref. for a comprehensive phenomenological review, [i^ . 



IV. SUMMARY 



The expansion parameter for heavy baryon xPT is TOtt/A^ as opposed to m^/A^ in the purely mesonic sector. 
Therefore it is crucial to calculate baryonic observables to at least NNLO to test the convergence of the heavy baryon 
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chiral expansion, as well as reduce the theoretical error bars of these calculations. To this order, we are forced to 
introduce a large number of LECs [see Ea. (|14ll5|) ]. which must be determined to have any predictive power, however, 
this is no small feat. To NNLO, there are more than ten LECs which must be fit to determine the nucleon masses [not 
including those associated with the A/A^ expansion, see Eq. and there are only two nucleons. Some of these 

can be fit from observables besides the mass, but even so, not all can be reliably determined. Lattice QCD simulations 
provide us with another means of determining the LECs, as we have the freedom to vary the quark masses, providing 
a much larger data set to fit these universal constants. Today, these lattice simulations are done with quark masses 
which are significantly larger than those in nature. It is therefore crucial to understand the quark mass dependence 
of hadronic observables: in particular for nuclear physics, those of the nucleons and deltas. In this way, the LECs 
can be fit, allowing the lattice data to be extrapolated from the valence and sea quark masses used to their physical 
values, thus making rigorous QCD predictions of these hadronic observables. 

In this work we have calculated the masses of the nucleons and the deltas in non-degenerate, two-flavor heavy 
baryon chiral perturbation theory to 0{rn^). One can obtain from the authors a Mathematica notebook containing 
the derived expressions in various limits, including the partially-quenched calculations worked out in Appendix 1X1 
relevant for current and future lattice data. We have highlighted the mass splittings of the nucleons and deltas as 
the 0{rfi?q^) contribution to these mass differences vanish, even away from the isospin limit. To LO these mass 
differences are proportional to (m^ — mu). Any deviations from this linear behavior are a signature of the NNLO 
mass contributions. This provides a clean handle to determine some of the LECs which contribute to the masses at 
O(m^), which can otherwise be rather difficult to extract. 

Another example of the utility of xPT is to predict the width of the deltas from the lattice. As the deltas are 
resonances, one cannot easily obtain their physical properties directly from real-valued Euclidean space simulations. 
As the lattice pion masses are brought down, one must deal with unstable particles. Instead, one can work with pion 
masses m^r > 300 MeV for which the deltas are stable and lattice calculations can be performed more simply than 
the unstable case (which also has yet to be explored). In this regime, one still hopes the chiral expansion is valid, ^ 
and uses the xPT expressions derived above, or the PQxPT expressions derived in Appendix IXI to fit the LECs and 
these consequently lead to predictions for the physical delta resonances. For example, their decay widths can be found 
to NNLO from knowing the pion mass and the constants /, A, A, um, cm, 7m, ^m, and qan- For simplicity the 
expression for the average delta width Qm(MT) in terms of these parameters is 



3m(MT) 



gAw / A 2 



0, 



r2 2 
' 24A 



[37m + 4(o-A/ - o-m) + 2a7\/] \, for < A 

for m^r > A 



(34) 



As a further application, the mass splitting between the nucleons and deltas can also be found from the various 
LECs. For the real part of the difference between the average nucleon and average delta masses, we have^ 

sRe(MT-Ms) = A-l-^tr(mg) -I-Btr(m^) +C[tr(TOg)]2 



I67r/2 

Aiv{mq) 

1 



Q 2 25 2 
'^9 A ~ ^.9aa 



^9 AN 



247r2/: 
3£(m^,/i) - 5gi^ 



2J^(m^, A,/i) 



5 



1 



1 



8Mf_ 



-E 



(4^/)2 



47r/ 



8Mb 



-G 



(35) 



Once the LECs are known from the calculation of baryon masses, say, one can use the chiral effective theory to calculate 
certain contributions to TriV — )■ ttN or ttA ttA scattering, for example, two processes which are complicated to 
explore on the lattice. 

The nucleons and deltas hold an important place in hadronic physics. As the lightest spin-i and spin-| baryons, 
they have received considerable attention theoretically to compare with a wealth of experimental data. Soon lattice 



^ This is only a conservative estimate as finite volume effects will modify the decay threshold of the deltas. 
Here we have used the following replacements A~G for particular combinations of the LECs: A = oim + '^(^M — f?Af ) + 7Af , B = + 

E=§g\ + 153i^ - f F = -l [fb^-^ + SbJ^ + + q^) + ^{tj + q^) + 3{t§ + q^)] , and G = fg\ + f - §gl^. 
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QCD calculations will enable an understanding of the nucleon and delta observables in terms of the quark masses. A 
crucial first step will be determining their masses from first principles. Current and foreseeable lattice calculations 
will rely upon chiral extrapolations. Our expressions for the masses in xPT and PQxPT enable such extrapolations, 
and provide a means to access physics that cannot be directly studied on current lattices. 
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APPENDIX A: PARTIALLY QUENCHED CHIRAL PERTURBATION THEORY 

The study of PQQCD is motivated by lattice QCD simulations. Presently all lattice QCD simulations are done 
with quark masses that are larger than those in nature, due to the computational cost of simulating lighter quarks. 
The valence quarks, those connected to external fields, can be treated independently from the dynamical sea-quarks 
that exist only in the closed quark loops. When the masses of the valence and sea quarks are taken to be different, 
one has partially quenched QCD (PQQCD), and in the limit that the valence and sea quark masses are equal, one 
has QCD. To exclude contributions from closed valence loops, one must include ghost quarks in the effective theory 
of the lattice action, which arc identical to the valence quarks, except that they have bosonic statistics. Thus the 
valence and ghost quark contributions to the PQQCD functional integral exactly cancel, leaving only the sea-quark 
contributions. 

As in QCD, one postulates that the vacuum of PQQCD spontaneously breaks the (graded) chiral symmetry down 
to the vector subgroup, and that one can perturb about the massless quark limit. The emerging low-energy effective 
theory of PQQCD is partially quenched xPT (PQxPT) [Ulillillillillil. Since PQQCD retains an axial anomaly, 
the flavor singlet field can be integrated out and there are no additional operators required that involve the singlet 
field. Moreover because the sea sector of PQQCD contains QCD, the same LECs of xPT appear in PQxPT with the 
same numerical values. Hence there is a means to determine xPT parameters using PQQCD lattice simulations. This 
hope has generated much activity in calculating baryon masses in partially quenched theories pll27ll50ll5lll5^ . 

In PQQCD, the quark part of the Lagrangian is 

6 

£ = E QM*^-"iQ)/Qfc- (Ai) 

This differs from the usual SU (2) Lagrangian of QCD by the inclusion of four extra quarks; two bosonic ghost quarks, 
(m, d), and two fermionic sea quarks, (j, I), in addition to the light physical quarks (w, d). The six quark fields transform 
in the fundamental representation of the graded SU{4:\2) group. They have been accommodated in the six-component 
vector 

Q'^ ^ {u,d,j,l,u,d). (A2) 

The quark fields obey the graded equal-time commutation relation 

(x)Qf (y) - (-l)''"'^Qf (y)Q?(x) = ^"^^..^^(x - y), (A3) 

where a, (3 and i,j are spin and flavor indices, respectively. Analogous graded equal-time commutation relations can 
be written for two Q's and two Q^'s. The grading factors 

take into account the different statistics of the quark fields of PQQCD. The quark mass matrix of 5?7(4|2) is given 

by 

mg = diag(m„,TO<i,TOj, m;,TO„,md). (A5) 
QCD is recovered in the limit rrij and m; — > m^. 
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1. Partially Quenched Lagrangian 

For massless quarks, the theory corresponding to the Lagrangian in Eq. (|A1|) has a graded SU{4\2)l <Xi 
SU{A\2)ji (g) U{l)v symmetry which is assumed to be spontaneously broken down to S'C/(4|2)y ® U{l)v in anal- 
ogy with QCD. The effective low-energy theory obtained by perturbing about the physical vacuum state of PQQCD 
is PQxPT 130,^^0,^^^. The emerging pseudo-Goldstone mesons have dynamics described at leading order 
in the chiral expansion by the Lagrangian 

£ = :^str (a^Sta^S) + A str (jngE^ + m^^s) (A6) 

where 

(A7) 

and the meson fields appear in 

$ = I I . (A8) 




The operation str() in Eq. ljA6p is a graded flavor trace. M and M are matrices containing bosonic mesons (with 
quantum numbers of qq pairs and qq pairs, respectively), while x ^^^d are matrices containing fermionic mesons 
(with quantum numbers of qq pairs and qq pairs, respectively). The upper 2x2 block of M contains the familiar 
pions and the remaining components are mesons formed with one or two sea quarks, see e.g. [5ll |. Expanding the 
Lagrangian Eq. (jA6|) to leading order, one finds that mesons with quark content QQ' are canonically normalized with 
their masses given by 

^QQ' = -^(^Q + ^Q')- ( A9) 

The flavor singlet fleld is defined to be $o = str(<i>)/\/2. As PQQCD has a strong axial anomaly U{1)a, the mass of 
the singlet field is on the order of the chiral symmetry breaking scale and has been integrated out of the theory [48l |. 
In this limit, however, the 77 two-point correlation functions deviate from their form in xPT. For a,b = u,d,j^l^u^d^ 
the leading-order rjaflb propagator with non-degenerate sea-quarks is 

^ iegSab i £a£b {q^ - m|^) {q^ ~ mfi) 

q^-ml+ie 2 {q^ - + le) {q^ ~ + le) {q^ ^ mj, + te) ' ^ ' 



where 



-1)1+"". (All) 



The mass irixy is the mass of a meson composed of (anti)-quarks of flavor x and y, while the mass mx is defined as 
= i (rrvjj + mfi'^. The flavor neutral propagator can be conveniently rewritten as 



Pa+eaebnab{Pa,Pb,Px), (A12) 



where 



Pn - ^ Pb = ^ P 



(7^ — mpi^ + le (7^ — m^^ + te q^ — mj^ 



Hab {A, B, C) 



-A- 



m 



X) 



B 



-C 



•'X 



(A13) 
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For our calculation, the nucleons and deltas must also be included in PQxPT '53,'5l'|. A straightforward way to do 
this is to use interpolating fields that have non-zero overlap with the nucleons and deltas. This leads to the tensor B^^'', 
which describes a 70-dimensional representation of SU(4|2), and the tensor T*^'^, which describes a 44-dimensional 
representation of SU(4|2), see [H^- The nucleon doublet ISfi is embedded in Bijk as 



ijk 



V6 



(AM) 



where the indices «, j, and k are restricted to 1 or 2. The deltas are contained simply as T'^^^ — T*-'*^ for i, j, and k 
restricted to 1 or 2. 

To write down the PQxPT Lagrangian in the baryon sector, we must also include the appropriate grading factors in 
the contraction of flavor indices. Such factors are included in the () notation defined in |5lll53l|. The free Lagrangian 
is given by 



£ = {Biv-D B) + 24^"^' {BBM) + 2(3['^'^' {BMB) 



(PQ) 



2(7ji^'3) ^gg) str(X) 



.f,j str(AI) . 

Above, the covariant derivative's action on either the B or the T fields is defined by 



(A15) 



(A16) 



and the operators A^, V^'' and Ai are all defined analogously to those in xPT. The free Lagrangian above contains 
more operators than the corresponding SU(2) Lagrangian. For this reason we have appended {PQ) superscripts to 
the coefficients. The relations between the free Lagrangian parameters of PQxPT and are found by matching. 
One finds 



CtM 

Cm 

7m 

CM 



3"m 



if:>{PQ) 

Q i-'m ' 



g"M 



1^' 



{PQ) 

M 



(A17) 

(A18) 

(A19) 
(A20) 



The Lagrangian describing the interactions of the 70 and 44 baryons with the pseudo-Goldstone particles is 

£ = 2a{BS^'BA^,) + 2p{BS''A^B) + 21^ (J" S^' A^,%) 

+ [(J^A^B) + (BA^T'')] . (A21) 



Matching this Lagrangian onto the SU(2) chiral Lagrangian in Eq. (|12|) . we find TL = gAA and C = 
more, we follow |5lj and choose to express a and /3 in terms of new parameters gA and gi, namely 

4 1 
a = -gA + gffi, 

R 2 1 



-.9A7V- 



Further- 

(A22) 
(A23) 



so that gA has the same value as in xPT. 

At higher orders in the chiral expansion, the situation in PQxPT parallels that of xPT in Sec. Recall that 
at higher orders, the Lagrangian can contain arbitrary functions of A/A^. We take this into account by implicitly 
treating the leading-order coefficients as functions of A/A^ expanded out to the required order. To maintain the 
Lorentz invariance of the theory, we use reparameterization invariance to generate the higher dimensional operators 
with fixed coefficients. In PQxPT the fixed coefficient Lagrangian is given by 



£ = - 



2Mb 



B^BvA 

. Mb 



'b^vAB 
. Mb 




Bv ■ A—- — B 



13 



+71 



(A24) 



As far as operators with unfixed coefficients, for the nuclcon mass we have the following partially quenched operators 
with two insertions of the axial-vector pion fields 

£ = _1_ jfof (^Q) [A ■ A): Bn,k + bt^'''^^ (^^yv.+v.Kri.+v.m'^^^ . B..,n 
^l^MPQ) Qk^j^^^^r^^ {A-A) + bl^^^^'^ W^' {v Av A)/' Bnjk 



^f^vAiPQ) (_)^,r,„+i^fe,.(^ . A)^"{v ■ A)}Bink + bf''^'^^ B'^^'B^j^Tt {vAvA)j, (A25) 
and the following operators with two insertions of the mass operator 

^5M(PQ) QkJ^l^^^^ ^^^^ ^ ^M(PQ) ■Qk3^^^ynl^^^^ ^j^^ 

^^M(PQ) str (X) + 6^ ^""^^ S'^^'^Byfe [str(Al)]2|. 

(A26) 

(A27) 
(A28) 
(A29) 
(A30) 

Carrying out the matching for the operators with coefficients b"^^^^^ and foj'^'-^'^' leads to the relations 



Matching the operators with coefficients 5*^^^'^'' above onto the SU(2) Lagrangian in Eq. 114|l . we find 



b{' = 


1 , Af (PQ) 


2 M(PQ) 


3 3 


l,Af(PQ) 
+ 2''4 




b¥ = 


2,A/(PQ) 


1 , Af(PQ) 


l,Af(PQ) 

6 3 + 


1 , A/(PQ) 
6'^4 + 


,M{PQ) 
"5 




l^M(PQ) _ 

2 3 


1 , Af(PQ) 


l,Af(PQ) 
36 + 


2,A/(PQ) 




bf = 


l^M(PQ) _ 

6 3 


6 4 + 


2,Af(PQ) 
36 + 


1 , A/(PQ) 

+ 


,M{PQ) 
"8 



uA,vA ^ ]_,A,vA(PQ) l,A^vA{PQ) _ \^A,vA{PQ) _^,A,vA(PQ) ,A,vA(PQ) . . ^ 

2 1 2 ^ 3^ 12" ^ ' ^ ' 



Higher-dimensional operators that contribute to the masses of the deltas in PQ^PT arise similarly from two 
insertions of the axial- vector pion fields 



r-tl 

'j'k 



+t^(^«^ (r^r^)str(A,A'') + if ^^^^rf (A^A.)f7:r,, + i^rf (A'').*'(A,)f7-r,,, 



+^3 - [T.T'^) str(A^A.) + r'^'iv-Av A^iT^f^, + if (v ■ A)f (^; • A)f 7-(;,, 



^A(PQ) 
,i;A(PQ) jTf 



(j^r^^) str(t; -AvA)), (A32) 
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and two insertions of the mass operator 



+4' {T^.T'') str(MM) + (T^MT'') str(M) + <f {T^X") [str(7WM)]2 i . (A33) 



The coefficients of these operators which do not have the superscript, (PQ), tf, i^, ig"^, and tf^ have the same 
numerical values as those in the SU(2) Lagrangian, Eq. p5|) . The operators with the (PQ) superscript are related to 
the SU (2) LECs by the following relations, 

,A,vA,A _ ^ ,A.vA,A(PQ) A,vA.A(PQ) r\oA\ 



Nucleons 



The relevant diagrams needed to obtain the nucleon mass to NNLO in PQxPT are depicted in |2y| and include 
hairpin contributions from the flavor-diagonal states. To leading order in PQxPT the nucleon masses are 



ivig — ^aj^i rrig 



1 



^ + 24^^^^ str(mQ), 



(A35) 



where the tree-level coefficients m'g and m'^ depend on the valence quark masses and are given in Table The 
next-to-leading contributions to the nucleon masses take the form 



M 



(3/2) _ 
B — 



29' 



AN 



(47r/) 



(A36) 



where T(m, S, /i) is given in Eq. H21(l . Additionally we have employed the abbreviations 



(A37) 



for contributions arising from the hairpin diagrams. The sums involving these two functions are over pairs of flavor- 
neutral states in the quark basis, e.g., above runs over "quVn, Vufld, and rjdrjd- In this way there is no double 
counting. The remaining sums in Eq. ljA36|l are over all loop mesons (j) of mass m^. The coefficients , A^^,, 
and B^^, appear in Table IVIl 

Finally the next-to-next-to-leading contributions to the nucleon mass are 



M 



(2) 



{Zb - l)Mi^> + 



A-Kf 



-6i ^'{m )b + -ba )b + -^h [mm )b + -64 ^'{mm Jg 



, ,M(PQ) , I 2n , l,Af(PQ) / , / X , 1,M(PQ) tl , , n , , M(PQ) , , , sVl 

-1-65 ^'str(mQ) -I- -6g ^'msStr(mQ) -I- -67 ^'m^ str(mQ) + 6g ^'[str(mQ)] 



2a 



(PQ) 



4,4,' 
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<t"p' 



+ 



i ST ( h^'MPQ) riB ,vA(PQ)j^B ,v 



4(4^/) 



£{■014,11) - ]-m% 



4(47r/)3 2-/ I, 1 -t^/3,</>0' + "2 -L^Q.^^'+Os + f)4 ii^^'+Og i^^^'j 



£{704,17141,^1) ~ -M {m4,m4') 



+ 



9 



8Ms (47r/) 



£{m4,fi) + ^m^ 



£{m4, 7041,11) + -A4*(m^,m0/) 



+ 



15 ffAJV 

8Ms (47r/)2 ■ 
eo-lf^'strfmc; 



£{704,11) + Y^n^,^. 



+ 



(47r/)2 

6gAJvo'Mstr(mQ) 
(47r/)2 



£{704,11) + -704 



44' 



- E ^^<p' 

44' 



£{014,014,, fi) + —M {704,014,) 



£{704,704, , ^) + -X^(m^,m0/) 



+ 



E-^? [^"("10, A, ^) + 70I] +^B^^, [J{ol4,7o^,,/:^,^l) + m'^ {704,70^,)] 



£{704, ^) + -704 



44' 



2 2 

£{704,704, ,11) + -M {704,704,) 



E '^f ^.P) + E <j'f<f>' \J{'m.4,704,,I^,[i) ^ JW'{t04,704,)\ 



(A38) 



The above expression is written quite compactly. It involves the tree-level coefficients m!^, m'g, {m'^)B, {m?)'g, 
{mm')B, and {771771')'^. These are listed in Table Ivl for the nucleon states labeled by B. The wavefunction renormal- 
ization Zb appearing in the above expression is given by 



Zb ~ 1 



(4^/)2 



E^f 



£{7714,^1) + ^7ol 



E 



2 2 

£{014,704,, fi) + -M {704,704,) 



3gAJV 

■ (4^/)2 



E [J{'m4,A,fi) + oil] + E [^('Ti^.m^/, A,/i) + A4^(m0,m0/)] 



(A39) 



Furthermore we have made use of abbreviations for the new non-analytic functions arising from loop integrals in 
Eqs. (jA38p and HA39|) . These are defined to be 



C{77i^,m^,,fi) = TL^,j,,[C{m^,^i),C{m^,,y),C{mx,y)], 
C{m^,m^',fi) = Ti4,4,,[C{m^,ix),C{m^,,^i),C{mx,tJ,)], 
J {771^, m^, ,5,fi)= 7^00' [J{m^ ,3,^), J{m^, ,5,fi), J{7nx , S, ^)] 



(A40) 



and arise from hairpin contributions. The various coefficients in the above sums over loop mesons and loop pairs of 
flavor- neutral mesons are listed in Tables IVlMXI One may check that in the limit TUj —> m„, m; rud the xPT 
results are obtained from the PQxPT expressions with the proper matching of coefficients, see Eq. (jA27|) - (jA34|l . 
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3. Deltas 

The relevant diagrams needed to obtain the delta masses at NNLO in PQxPT are depicted in fi^. These include 
hairpin contributions from the flavor-diagonal propagator. To leading order in PQ^PT the delta masses are 

= g7M mr - 2 cta/ str(mQ), (A41) 



where the coefficients ttit are listed for deltas T in Table ITU At next-to-leading order, we have 



M. 



(3/2) 



72^/2 



9\n 

(4^/)2 



4,4,' 



(A42) 



where ^-"(771, A, /i) is given by Eq. (|21|l and the functions A^"(m0,TO0/) and ^{m^^m^i ^5, ^) appear in Eq. ljA37p . 
The coefficients A^, A^^/, -Bj and Sj^, appear in Table IXll 

At next-to-next-to leading order, we have the following contribution to the delta masses 



M, 



+ 



+ 



1 

47^7 

2 7M 
(47r/)2 

(4^ 



-^{m^)T + -^{mm')T + str{m%) + -^rriT str(mQ) + ig^ [str(mQ)]^ 



3 ' ' 3 

0<^' 



4(47r/) 



4(4^/) 



5 ffAA /(T' 



15 ffAJV 

16 {'iixfYMB 



4,4,' 



lOgAAO'Mstr(mQ) 
3(47r/)2 



E^ 



E^^ 



3flAJv'^M'^'str(mQ) 
(47r/)2 



44' 



+ 



lOgAAlM 

9(47r/)2 



E^; 



T 
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3ffAiV 

'2(4^/)2 



J(m0, - A, /i) + ^ Gij^0/ J(m0, m^, , -A, 



(A43) 



As with the NNLO contributions to the nucleon masses, the above expression is written quite compactly. It involves 
the tree-level coefRcients niT, (w^)t, and {mm')T which are listed in Table ITU for the delta states T. The NNLO 
result involves the wavefunction renormalization Zt, which we find to be 



Zt - 1 



5gAA 

'3(47r/)2 



3flAJV 

"2(47r/)2 



4,4,' 



26 



^ BlJ{m4, - A, m) + E -5^0' J(»n0, m^/ , - A, ^) 



(A44) 



The functions arising from hairpin contributions are listed in Eq. (|A40|I . The various coefficients in the above sums 
are listed in Tables IXH IXVI Finally, one can check that the xPT result is recovered in the limit nij rriu, mi rrid, 
and by utilizing the matching relations given in Eqs. IX27ll . (TOTTl . IMall and ICKM . 
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TABLE V: The tree-level coefficients in xPT and PQxPT. The coefficients ms, m's, m'g, (m^)s, (m^)'B, {mm')B, and {mm')'B 
are listed for the nucloons B. 





TUB 


niB 




(m2)B 




{mm')B 


{mm'YB 


p 


ruu 


rUu + 2md 


5mu + rud 




5mu + ml 


ml + 5mumd 


ml — Amumd 
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Sm^md + ml 


-Am^md + ml 
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APPENDIX B: PARTIALLY QUENCHED COEFFICIENT TABLES 

Here we list the tables needed for the partially quenched calculation which have not already been listed in the text. 
One can obtain from the authors a Mathematica notebook which evaluates PQxPT and xPT expressions for the 

nucleon and delta masses in various limits. Wc consider the case of non-dcgcncratc valence and sea quarks, as well as 
non-degenerate valence and degenerate sea quarks and further degenerate valence and sea quarks. 
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TABLE VI: The coefficients Af, Af^,, Bf, and Sj^, in PQxPT. Coefficients are listed for the nucleons B, and for Af and 
Bf are grouped into contributions from loop mesons with mass m<^, while for A'^^, and B^^, are grouped into contributions 
from pairs of quark-basis r)q mesons. 









/I? 














p 


UdA + 25A51 + flf /4) - ffAfli - 5flf /4) 




n 





KflA - ffA5i - 55?/4) + 2ffAffi + S?/4) 






ju 




jd 


Id 


P 


3[^9a + 9a9i + 9iI^) 


3(23i+flAfli+fl?/2) fl?/4 


3f/4 


n 




ff?/4 


\{2g\ + gA9i + 9f/'2) U^gl + 9a9i + ff?/2) 














iji, iji, 


'/„ '/J 


11,111,, 




P 


9a + 9A91 + 5i/4 


SA5i + c/?/2 






n 


S?/4 


5Affi + 5? /2 


9a + 9A91 + ff?/4 
















rtu 




Vd 




P 


1 

9 


6 
9 







n 





6 
9 


1 

9 






ju 




jd 


Id 


P 


1 

9 


1 
9 


2 
9 


2 
9 


n 


2 
9 


2 
9 


1 
9 


1 
9 


















VdVd 




P 


2 

!) 


4 
9 


2 
9 




n 


2 
9 


4 
9 


2 
9 





20 



TABLE VII: The coefficients , C, 
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c. 



Dtf,, and D^^i in PQxPT. The coefficients C^, 
and D^^i arc identical for all nucloons and deltas. The remaining coefficients are listed for nucleon states B, and 



for C^, and 



are grouped into contributions from loop mesons with mass m,/,, while for C^^, and D^^, are grouped into 
contributions from pairs of quark-basis rjq mesons. If a particular meson or pair of flavor-neutral mesons is not listed, then the 
value of the coefficient is zero for the nucleons. 
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TABLE VIII: The coefficients , Ef^,, E'^ and E'^^, in PQxPT. Coefficients are Usted for the nucleons, and for E^ and E'^ 



are grouped into contributions from loop mesons with mass m^, while for E^^i and E^ 
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TABLE IX: The coefficients and F^^, in PQxPT. Coefficients are listed for the nucleons, and for Fj, are grouped into 
contributions from loop mesons with mass m^, while for F^^, are grouped into contributions from pairs of quark-basis r}q 
mesons. 
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TABLE X: The coefficients Gf and Gf^, in PQxPT. Coefficients are listed for the nucleons, and for are grouped into 
contributions from loop mesons with mass m^, while for G^^/ are grouped into contributions from pairs of quark-basis rjq 
mesons. 
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TABLE XL The coefficients 4 J, A^^,, and Sj^, in PQxPT. Coefficients are listed for the delta states T, and for and 
are grouped into contributions from loop mesons with mass mc/,, while for A'^^, and B^^, are grouped into contributions 
from pairs of quark-basis rjq mesons. 
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TABLE XII: The coefficients Gj, Gj^,, and Dj^, in PQxPT are listed for the deha states T, and for Gj and L>J are 
grouped into contributions from loop mesons with mass m^, while for Gj^/ and D"^^, are grouped into contributions from pairs 
of quark-basis riq mesons. If a particular meson or pair of flavor-neutral mesons is not listed, then the value of the coefficient is 
zero for all deltas. The coefficients C4,, C 4,4,1, D4, and D^^i axe identical for all nucleons and deltas and appear in Table IVTTl 
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TABLE XIII: The coefficients FJ and FJ'^, in PQxPT. Coefficients are listed for the delta states T, and for Fj, are grouped 
into contributions from loop mesons with mass m^, while for F^^, are grouped into contributions from pairs of quark-basis riq 
mesons. 
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TABLE XIV: The coefficients and E^^, in PQxPT. Coefficients are listed for the delta states T, and for F^, are grouped 
into contributions from loop mesons with mass m^, while for FJ^/ are grouped into contributions from pairs of quark-basis m 
mesons. 
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TABLE XV: The coefficients and Gj^, in PQxPT. Coefficients are listed for the delta states T, and for G^ are grouped 
into contributions from loop mesons with mass m^, while for Gj^/ are grouped into contributions from pairs of quark-basis r/q 
mesons. 
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